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Abstract 

Background: GATA transcription factors are essential for self-renewal of the small intestinal epitheliunn. Gata4 is 
expressed in the proxinnal 85% of small intestine while Gata6 is expressed throughout the length of small intestine. 
Deletion of intestinal Gata4 and Gata6 results in an altered proliferation/differentiation phenotype, and an up-regulation 
of SAM pointed domain containing ETS transcription factor (Spdef), a transcription factor recently shown to act as a 
tumor suppressor. The goal of this study is to determine to what extent SPDEF mediates the downstream functions of 
GATA4/GATA6 in the small intestine. The hypothesis to be tested is that intestinal GATA4/GATA6 functions through 
SPDEF by repressing Spdef gene expression. To test this hypothesis, we defined the functions most likely regulated by the 
overlapping GATA6/SPDEF target gene set in mouse intestine, delineated the relationship between GATA6 chromatin 
occupancy and Spdef gene regulation in Caco-2 cells, and determined the extent to which prevention of Spdef 
up-regulation by Spdef knockout rescues the GATA6 phenotype in conditional Gata6 knockout mouse ileum. 

Results: Using publicly available profiling data, we found that 83% of GATA6-regulated genes are also regulated by 
SPDEF, and that proliferation/cancer is the function most likely to be modulated by this overlapping gene set. In 
human Caco-2 cells, GATA6 knockdown results in an up-regulation of Spdef gene expression, modeling our mouse 
Goto6 knockout data. GATA6 occupies a genetic locus located 40 kb upstream of the Spdef transcription start site, 
consistent with direct regulation of Spdef gene expression by GATA6. Prevention of Spdef up-regulation in conditional 
Gata6 knockout mouse ileum by the additional deletion of Spdef rescued the crypt cell proliferation defect, but had 
little effect on altered lineage differentiation or absorptive enterocytes gene expression. 

Conclusion: SPDEF is a key, immediate downstream effecter of the crypt cell proliferation function of GATA4/GATA6 in 
the small intestine. 
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Background 

The mature mammalian small intestine is lined by a highly 
specialized epithelium that regenerates itself in a tightly 
controlled manner resulting in a lineage distribution and 
gene expression patterning that is perfectly suited for the 
absorption of nutrients. The epithelium is organized into 
crypt-villus structures in which the Crypts of Lieberkiihn 
contain stem cells that produce proliferating, transit- 



* Correspondence: stephen.krasinski@childrens.harvard.edu 
^Division of Gastroenterology and Nutrition, Department of Medicine, 
Children's Hospital Boston, and Harvard Medical School, 300 Longwood 
Avenue, Boston, MA 021 15, USA 

Full list of author information is available at the end of the article 

(3 BioMed Central 



amplifying (TA) cells that differentiate into five principal 
post-mitotic cell types comprised of one type of absorptive 
cell (absorptive enterocytes) and four types of secretory 
cells (enteroendocrine, goblet, Paneth, and tuft cells [1]). 
Absorptive enterocytes express digestive enzymes and 
transporters necessary for the terminal digestion and 
absorption of nutrients. Mucus-secreting goblet cells 
and defensin-secreting Paneth cells maintain a dynamic 
mucosal defensive barrier. Enteroendocrine cells secrete 
hormones that regulate gastrointestinal processes. Tuft 
cells, recently shown to be an independent secretory lineage 
[2], secrete opioids and produce enzymes that synthesize 
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prostaglandins, suggesting a role in inflammation. Absorp- 
tive enterocytes, goblet cells, enteroendocrine cells, and 
tufl: cells migrate up the crypt to populate the villi, whereas 
Paneth cells migrate to the base of crypts. The differenti- 
ated cells eventually undergo apoptosis and are shed into 
the lumen. Cells of the villus epithelium turn over in 
3-4 days, whereas Paneth cells at the base of crypts turn 
over at a slower rate of 3-6 weeks. 

Current models of intestinal epithelial renewal suggest 
that long-lived, multipotent stem cells produce progenitors 
that undergo a series of transitions that ultimately give rise 
to the individual cell lineages [1,3]. The Wnt, hedgehog, 
and bone morphogenetic protein signaling pathways regu- 
late intestinal proliferation and differentiation, while the 
Notch signaling pathway plays a central role in determin- 
ing epithelial cell fate. The first decision selects absorptive 
vs. secretory progenitors. Activated Notch signaling re- 
sults in the transcriptional activation of its principal 
intestinal target, hairy and enhancer of split 1 {Hesl), 
which encodes a transcription factor that selects the 
absorptive enterocyte lineage. Progenitor cells that escape 
Notch signaling and activation of Hesl gene transcription 
express atonal homolog 1 {Atohl, formerly called Mathl), 
which encodes a transcription factor that selects the 
secretory cells. Additional regulators that function in 
secretory cell differentiation include: growth factor in- 
dependent 1 {Gfil) that distinguishes enteroendocrine 
from goblet/Paneth progenitors; neurogenin 3 {NeurogS) 
that specifies the enteroendocrine lineage; SAM pointed 
domain-containing Ets transcription factor {Spdef), a GFIl 
target, that promotes goblet differentiation; and SRY-box 
containing gene 9 {Sox9) and ephrin type B receptor 3 
(EphbS), Wnt targets that are necessary for the differen- 
tiation of Paneth ceUs and their localization to the crypt 
base, respectively. 

Recently, we showed that members of the GATA family, 
an ancient family of transcription factors that bind WGA- 
TAR motifs in DNA, play essential roles in crypt cell 
proliferation, secretory cell differentiation, and absorptive 
enterocyte gene expression. Gata4 and Gata6 are expressed 
in the intestinal epithelium, but whereas Gata6 is expressed 
throughout the length of the small intestine, Gata4 is 
expressed in the proximal 85% of small intestine and is 
sharply down-regulated in the distal ileum [4-6]. Using 
conditional knockout technology, we [5,7] and others [8] 
have shown that GATA4 functions to promote a 'jejunal' 
pattern of absorptive enterocyte gene expression and func- 
tion while repressing an lleaF pattern. Using single and 
double conditional knockout approaches for Gata4 and 
Gata6y we found that in the ileum of single Gata6 con- 
ditional knockout mice, where Gata4 is not normally 
expressed, or throughout the small intestine of double 
Gata4/Gata6 conditional knockout mice, crypt cell prolif- 
eration and enteroendocrine cell specification are decreased. 



Paneth cells are replaced by a goblet-like ceU type, and 
the expression of specific absorptive enterocyte genes is 
altered [4]. We also noted that Spdef^ a transcription factor 
expressed in secretory progenitors, goblet cells and Paneth 
cells that functions in goblet and Paneth cell differentiation 
[9,10], was up-regulated [4]. Using a conditional Spdef 
over-expression model, Noah et al. [10] described an intes- 
tinal phenotype that, with the exception of the changes in 
absorptive enterocyte gene expression, essentially pheno- 
copies that of our Gata4/Gata6 conditional knockout mice: 
crypt cell proliferation is decreased, enteroendocrine and 
Paneth cells are decreased, and goblet-like cells are in- 
creased. GATA6 is co-expressed with SPDEF in the same 
lineages in the small intestine [4]. Based on these findings 
[4], we hypothesized that GATA4/GATA6 regulates crypt 
cell proliferation and secretory cell differentiation in the 
small intestine by repressing Spdef gene expression. To 
test this hypothesis, we defined the functions most likely 
regulated by the overlapping GATA6/ SPDEF target gene 
set in mouse intestine, delineated the relationship between 
GATA6 chromatin occupancy and Spdef gene regulation 
in Caco-2 cells, and determined the extent to which 
prevention of Spdef up-regulation by Spdef knockout 
rescues the GATA6 phenotype in conditional Gata6 knock- 
out mouse ileum. 

Results and discussion 

GATA6 and SPDEF regulate similar subsets of genes 

To gain insight on the relationship between GATA6 and 
SPDEF in the small intestine, we scanned the overlap of 
gene targets using publicly available gene profiling data 
from conditional Gata6 and Spdef knockoMt mouse intes- 
tine [4,9]. Previously, we identified 2564 genes whose 
expression is altered in ileum by conditional Gata6 de- 
letion [4]. Network analysis of this gene set indicated 
an up-regulation of p53 targets and a down-regulation of 
c-MYC targets (Additional file 1: Figure SI), consistent 
with a decrease in cellular proliferation. Of the 2564 genes 
altered by conditional GataG knockout, 83% (2119) were 
also altered by Spdef knockout (Figure 1 A), a far greater 
overlap than would be expected from a similar-sized, ran- 
domized allocation of genes (P < 10'^^, Fishers Exact Test). 
The changes in expression of this subset when Gata6 is 
deleted were analyzed by gene set array analysis. Using 
Database for Annotation, VisuaUzation and Integrated 
Discovery (DAVID), we conducted functional annotation 
clustering of all the major pathways (listed in the KEGG, 
Biocarta and BBID databases). We found that Wnt signal- 
ing was the function most likely to be effected (Figure IB), 
consistent with regulation of crypt cell proliferation. Using 
Gene Set Enrichment Analysis (GSEA), a publically 
available bioinformatics tool that delineates gene expres- 
sion data for enrichment of pre-defined gene-sets, APC 
target network was one of the top three networks affected 
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Figure 1 Overlap of genes in small intestine modulated by conditional Gata6 or Spc/ef knockout. (A) Publicly available gene profiling data 
from conditional Gata6 (GSE22416) and Spdef (GSE14892) knockout mouse intestine show a significant overlap. (B) Functional annotation 
clustering of pathways shows that the top-enriched cluster, in overlapping Gota6 and Spdef genes, contains the Wnt pathway as its main function 
(ES score 2.58). (C) GSEA analysis shows, in the overlapping Goto6 and Spdefgene segment, APC targets as the second most highly enriched 
cluster (ES score 1.57). 



by Gata6 deletion with a Normalized Enrichment 
Score (NES) of -1.57 (Figure IC). This analysis reveals a 
very strong overlap in gene targets between GATA6 and 
SPDEF, and suggests that the principal function of this 
overlap involves cellular proliferation. 

GATA6 regulates Spc/ef gene expression and occupies a 
locus in the Spdef 5'-flanking region in human Caco-2 cells 

Conditional Gata6 deletion in mice produces an up- 
regulation of Spdef g^ne expression in ileum while knock- 
out or over-expression of Spdef has no effect on Gata6 
expression [4], suggesting that Spdef is regulated down- 
stream by GATA6. To determine the extent to which 
this process is conserved, and to further explore the role 
of GATA6 in regulating Spdef gene expression, we char- 
acterized the effect of Gata6 knock-down on Spdef gene 
expression in Caco-2 cells. The Caco-2 cell line is a human 
colorectal adenocarcinoma-derived cell line that expresses 
GATA6, but very littie GATA4 [11], similar to the ileum. 



We screened five GATA6 short-hairpin RNA (shRNA) 
knockdown lentiviral constructs, and found that only one 
resulted in a statistically significant knockdown of Gata6 
mRNA (55%, P < 0.05, Figure 2A) and a concomitant 
decrease in GATA6 protein (Figure 2B). Spdef mRNA 
was up-regulated 2-fold in Gata6 knock-down {G6kd) 
cells (P < 0.05, Figure 2A), demonstrating that the Caco-2 
cells model the up-regulation of Spdef gene expression 
observed in conditional Gata6 knockout mouse ileum [4]. 

To gain insight on the relationship between GATA6 and 
Spdef gene expression, we examined GATA6 chromatin 
occupancy at the Spdef gene locus. Because we were 
unable to immunoprecipitate bound chromatin from 
crosslinked mouse intestinal epithelial cells using existing 
antisera, we utilized the human Caco-2 intestinal cell 
culture model in which GATA6 ChIP assays have been 
performed previously [12]. Using a publicly available 
GATA6 chromatin immunoprecipitation-high throughput 
sequencing (ChlP-seq) database in Caco-2 cells [12], we 
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Figure 2 GATA6 regulates Spcfef expression and occupies an enhancer in the Spc/ef 5'-flanl<ing region in Caco-2 cells. (A) Quantitative 
RT-PCR analysis sliowing reduced expression of Goto6 mRNA, and enlianced expression of Spdef roHNA in Caco-2 cells infected with a lentivirus 
vector expressing an shRNA for human Goto6 mRNA (G6/cd) (mean ± SD, n = 5, ""P < 0.05). An shRNA vector for GFP was used as a control (Ct/). 
(B) Western analysis showing reduced abundance of GATA6 in G6kd Caco-2 cells. B-actin was used as an internal loading control. (C) Schematic 
representation of the human Spdef 5'-flan king region showing GATA6 occupancy at a locus -40 kb upstream of the transcription start site (TSS). 
GATA6 sequence tag density is shown as a 'wiggle file', and a statistically significant GATA6-occupied locus was defined by MACS peak analysis 
[29] {dotted box). (D) ChIP assays on chromatin obtained from Caco-2 cells using a GATA6 antibody and three sets of overlapping primers centered on 
the GATA motif showing increased GATA6 occupancy (mean ± SD, n = 4, *P < 0.05). The salivary amylase-ala {Amy-l) TSS was used as a negative 
control. Chip assays using anti-rabbit IgG was used as a control for non-specific immunoprecipitation and primer efficiency. 
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identified a GATA6-occupied c/^-regulatory region ap- 
proximately 40 kb upstream of the Spdef transcription 
start site (TSS) that mapped to the Spdef gene based on 
closest distance to any TSS (Figure 2C). To confirm and 
localize GATA6 occupancy at this locus, we performed 
GATA6 Chip assays using a multiple, overlapping primer 
strategy with amplicons of -200 bp per primer pair 
(Figure 2C) centered at the only GATA motif within 
this region. We found increases of all three amplicons with 
amplicon from primer set 2 being significantly greater 
(-70%, P < 0.05, n = 4) than the Amyl TSS negative control 
(Figure 2D). Although our ChIP data reveal a modest in- 
crease in enrichment, it is nonetheless statistically signifi- 
cant and thus confirmatory of ChlP-seq data conducted 
by others [12] (see MACS peak data. Figure 2C). We 
also conducted an IgG control ChIP assay to show that 
enrichment with primer set 2 was not due to differences 
in primer efficiency from that of the Amyl TSS, or to 
non-specific antibody binding (Figure 2D). These data 
confirm previous ChlP-seq analysis that GATA6 occupies 
a c/5-regulatory region that maps to the Spdef gene [12], 
and, together with Spdef up-regulation in G6kd cells, is 
consistent with the notion that GATA6 directly represses 
Spdef gene transcription. While mapping of occupancy 
sites to the nearest TSS is a well recognized method for 
defining transcription factor targets on a global basis 
[12,13], direct regulation of Spdef gene expression by 
GATA6 at this site will need to be confirmed using chromo- 
some conformation capture or other techniques [14]. 

Transcriptional repression is highly complex, especially in 
eukaryotes, generally involving the recruitment of specific 
co-repressors and the local modification of histone tails and 
chromatin structure [15]. GATA factors are well known to 
mediate gene repression, and the mechanisms are begin- 
ning to be understood. GATA4 has been shown to directly 
repress cardiac genes [16] while GATAl has been shown 
to directly repress hematopoietic genes [17,18]. Both 
GATA4 and GATAl interact directly with corepressor 
complexes including the nucleosome remodeling and his- 
tone deacetylase (NuRD) complex and the polycomb 
repressive complex 2 (PRC2) [16], and are necessary 
for the multiple modifications of histone tails, but the 
extent to which GATA factors recruit co-repressors and/ 
or modulate histone tails in the small intestine remains to 
be determined. 

Spdef knockout rescues the proliferation defect in 
conditional Gata6 knockout mice 

We next asked to what extent Spdef deletion rescues the 
conditional Gata6 knockout phenotype by analyzing single 
and double GataG/Spdef knockout mice. In previous 
studies, conditional deletion of Gata6 [4] or deletion of 
Spdef [9] resulted in greatly diminished levels of Gata6 or 
Spdef m^^K, respectively, in the intestine that correlated 



with reduced protein levels and altered intestinal pheno- 
types. Analysis of mRNA in the present study showed that 
Gata6 and Spdef As were both expressed in Ctl mice, 
Gata6 mRNA was nearly undetectable while Spdef mR^ A 
was expressed in Gata6AIE mice, Gata6 mRNA was 
expressed while Spdef mRNA was nearly undetectable 
in SpdefKO mice, and both GataG and Spdef vc\!^ks> were 
nearly undetectable in the double knockout {DKO) mice 
(Additional file 1: Figure S2), verifying our models. 

Crypt cell proliferation is essential for the continuation 
of intestinal epithelial renewal. One of the principal pheno- 
typic outcomes of conditional Gata4/Gata6 deletion is 
a -30% reduction in the number of Ki67- and BrdU-positive 
cells, and a concomitant decrease in villus height and vil- 
lus epithelial cell number [4] resulting in a reduction in 
absorptive surface area. To define the role of SPDEF in 
mediating the decrease in crypt cell proliferation when 
Gata6 is conditionally deleted, we stained ileal segments 
for Ki67 and BrdU (Figure 3A), markers of the non-Go 
and S -phases of the cell cycle, respectively. The number of 
Ki67-positive and BrdU-positive crypt cells was signifi- 
cantly reduced -40% in Gata6AIE mice (Figure 3B), as 
previously reported [4], but unchanged in SpdefKO mice, 
as compared to controls. The number of Ki67- and 
BrdU-positive crypt cells in DKO mice was also similar 
to controls, indicating that the additional deletion of 
Spdef rescues the decrease in crypt cell proliferation 
observed in Gata6AIE mice. These data support the notion 
that GATA6 maintains crypt cell proliferation by down- 
regulating Spdef gene expression. 

SPDEF does not regulate the secretory cell differentiation 
function of GATA6 

Previously, conditional deletion of Gata6 in the ileum 
resulted in a 25-40% reduction in the number of chromo- 
granin A (CHG A) -positive cells, and in the mRNA abun- 
dances of Chga and NeurogS [4], consistent with a decrease 
in enteroendocrine lineage commitment. Though not sta- 
tistically significant, we found a 15-50% reduction in the 
number of CHGA-positive cells (Figure 4A), and in the 
mRNA abundances of Chga and NeurogS (Figure 4B) in the 
Gata6AIE mice, in general agreement with our previous 
study [4]. We further found that although the pattern for 
these three measurements in the SpdefKO mice was similar 
to that in Ctl mice, consistent with previous data for distal 
intestine [9], the pattern in the DKO mice was similar to 
that in the Gata6AIE mice (Figure 4A and B). Together, 
these data show that enteroendocrine cell commitment 
is reduced in both Gata6AIE and DKO mice, indicating 
that GATA6 promotes enteroendocrine cell commitment 
independently of SPDEF. 

Previously, we also showed that conditional Gata6 
deletion resulted in a transformation of Paneth cells into 
Mucin-2 (MUC2) -enriched goblet-like cells at the base of 
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Figure 3 Crypt cell proliferation is decreased in Gata6AIE mice. (A) Immunostaining in ileum for the proliferation markers Ki67 {top row) and 
BrdU {bottom row). (B) Quantification of Ki67- and BrdU-positive cells/crypt in each group. Cells were counted as described in Methods. Gata6AlE 
mice had significantly fewer Ki67- or BrdU-positive cells than each of the other three groups (P < 0.01 in each case), consistent with a decrease in 
crypt cell proliferation. 



crypts [4]. In the present study, immunostaining for 
CRS4C, a Paneth-specific marker, was slightly reduced 
in SpdefKO mice as compared to controls, consistent 
with previous data [9], but was greatly reduced in both 
Gata6AIE and DKO mice as compared to controls 
(Figure 5A). The number of crypt cross-sections with 
at least one CRS4C-positive cell was nearly 100% in Ctl 
and SpdefKO mice, but was less than 30% in both 
Gata6AIE and DKO mice (P < 0.01 for each. Figure 5B). 



The mRNA abundance for lysozyme {Lyz), another Paneth 
marker, was not significantly different between Ctl and 
SpdefKO mice, but was significantly reduced in Gata6AIE 
and DKO mice as compared to controls (P < 0.01 for each, 
Figure 5C). PAS staining was greatly increased in crypts 
of Gata6AIE mice, as previously shown [4], noticeably 
reduced in SpdefKO mice, especially in villus goblet 
cells, as previously shown [9], and similar to controls in 
DKO mice (Figure 5A). MUC2 immunostaining was greatly 



CHGA-positive cells Chga nnRNA abundance Neurog nnRNA abundance 




Figure 4 Enteroendocrine cell allocation is decreased in Gata6AIE and DKO mice. (A) Quantification of CHGA-positive cells in each group. 
Sections of ileum were stained for CHGA, and the number of positive cells/1000 epithelial cells was determined as described in Methods. (B) Chgo 
and Neurog3 mRNA abundance in each group. Qata6AIE and DKO mice had generally lower numbers of CHGA-positive cells and Chga and 
NeurogS mRNA abundances than Ctl mice, consistent with a decrease in enteroendocrine lineage allocation. 
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Figure 5 Paneth cells are decreased In Gata6AIE and DKO mice. (A) Immunostaining in ileum for CRS4C {first row) and MUC2 {third row), and 
cliemical staining using tine PAS reaction {second row). (B) Quantification of CRS4C-positive cells in each group. Sections of ileum were stained for 
CRS4C, and the number of positive cells/crypt cross section was determined as described in Methods. (C) Lyz and Muc2 mRNA abundance in each 
group. GotoGAiE and DKO mice had significantly lower numbers of CRS4C-positive cells and Lyz mRNA abundance than Cti mice, consistent with 
a decrease in mature Paneth cells. 



enriched in crypts of Gata6AIE mice, consistent with our 
previous report [4], but similar to controls in SpdefKO and 
DKO mice (Figure 5A). Generally consistent with this 
observation, Muc2 mRNA abundance in Gata6AIE mice 
was 2.5-fold higher as compared to controls, though not 
statistically significant, whereas Muc2 mRNA abundances 
in SpdefKO and DKO mice were similar to controls 
(Figure 5C). These data indicate that the additional dele- 
tion of Spdef did not rescue the decrease in the terminal 
differentiation of Paneth cells observed in Gata6AIE 
mice, but did at least partially prevent their conversion 
to MUC2-enriched goblet-like cells. Thus, SPDEF is not 



necessary for the GATA-mediated differentiation of Paneth 
cells, but does function in their default differentiation into 
goblet-like cells in the absence of Gata6, This is consistent 
with the general function of SPDEF to promote the differ- 
entiation of goblet cells [9,10]. 

In spite of the common secretory cell phenotypes in 
Gata6 deletion and over-expression models [4,10], 

our data show that Spdef delQiion did not rescue the Gata6 
knockout defects in enteroendocrine lineage commitment 
or in Paneth cell differentiation. Closer scrutiny suggests 
that the secretory cell phenotypes in the Gata6 deletion 
and over-expression models are not identical. While 
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the decline in enteroendocrine and Paneth cells, and accu- 
mulation of goblet cells, appear similar, the underlying 
alteration in Paneth and goblet cell differentiation is 
different. Conditional Gata6 deletion has no effect on 
the normal commitment and differentiation of goblet 
cells on villi but reveals a conversion of Paneth cells at 
the base of crypts into a goblet-like cell type. These cells 
do not express defensins, but express abundant Muc2; 
they also express abundant Sox9 and EphbS [4], indicating 
that they are committed and targeted Paneth cells. On the 
other hand, conditional over-expression results in a 

generalized increase in goblet cells at the expense of all 
other epithelial cell types [10]. Hence, these mice show a 
decline in Paneth cell specification rather than a defect in 
the terminal differentiation of committed Paneth cells, as 
observed in our conditional Gata6 knockout model. 

SPDEF does not regulate the absorptive enterocyte gene 
expression function of GATA6 

We next examined whether Spdef loss affected GATA6- 
dependent absorptive enterocyte gene expression. Condi- 
tional Gata6 deletion resulted in a down-regulation of 
specific absorptive enterocyte genes in ileum that in- 
clude lipid transporters and apolipoproteins, and an up- 
regulation of genes in absorptive enterocytes normally not 
expressed or expressed at low levels in small intestine, but 
expressed at high levels in colon [4]. Using marker genes 
for these two patterns, apolipoprotein Al {Apoal) and 
carbonic anyhdrase 1 {Carl), respectively, we found that 
their mRNA abundances were down-regulated and up- 
regulated, respectively, in Gata6AIE mice (Figure 6), 
consistent with our previous findings [4]. Apoal and 
Carl mRNA abundance in SpdefKO mice was similar 
to controls, whereas that in DKO mice was similar to 
GataSAIE, indicating that SPDEF does not play a role 



in mediating the conditional Gata6 knockout alteration in 
absorptive enterocyte gene expression. 

GATA6 is expressed in crypts and in mature absorptive 
enterocytes [4,5], whereas SPDEF is expressed in secretory 
progenitors, and in goblet and Paneth cells, but is not 
expressed in mature absorptive enterocytes [9,10]. While 
it is possible that SPDEF could instruct gene expression in 
absorptive enterocytes through a process that originates in 
progenitors early in the differentiation process, we found 
that SPDEF did not regulate the GATA6 targets studied 
here. As in the GATA4-specific pathway, in which GATA4 
regulates its targets in mature absorptive enterocytes on 
villi rather than in crypt progenitors [19], we believe that 
GATA6 also regulates absorptive enterocyte gene expres- 
sion within mature absorptive enterocytes on villi. 

Conclusion 

Previously, we defined two fundamental pathways of GATA 
regulation in the small intestine, one mediated exclusively 
by GATA4 (GATA4-specific pathway), and one regulated 
by GATA4 or GATA6 (GATA4/GATA6-redundant path- 
way) (Figure 7). In the GATA4-specific pathway, GATA4, 
but not GATA6, activates and represses a subset of absorp- 
tive enterocyte genes, and by virtue of its expression in the 
proximal 85% of small intestine and lack of expression in 
distal ileum [7], distinguishes proximal intestinal from distal 
ileal gene expression and function [5,7,8]. For the GATA4/ 
GATA6-redundant pathway, either GATA4 or GATA6 
(which is expressed throughout the small intestine, includ- 
ing distal ileum) promote intestinal epithelial renewal by 
supporting crypt cell proliferation, enteroendocrine lineage 
commitment, Paneth cell differentiation, and absorptive 
enterocyte gene expression [4]. Here, we show that the 
maintenance of crypt cell proliferation function of the 
intestinal GATA4/GATA6-redundant pathway is dependent 
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Figure 6 Absorptive enterocyte gene expression is altered in Goto6AIE and DKO mice. Goto6AIE and DKO mice had significantly lower 
Apoal, and significantly higher Corl mRNA abundances than Ctl mice. 
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Figure 7 Model for the known GATA pathways in mature small 
intestine, and the placement of SPDEF within those pathways. 

GATA4 (purple), but not GATA6 [blue], activates and represses a subset of 
absorptive enterocyte genes (GATA4-specific), wliereas eitlier GATA4 or 
GATA6 ipurple/blue) promote crypt cell proliferation, enteroendocrine 
lineage commitment, Paneth cell differentiation, and absorptive 
enterocyte gene expression (GATA4/GATA6-redundant). GATA4/GATA6 
promotes the crypt cell proliferation function by repressing the gene 
encoding SPDEF (green). 



on SPDEF (Figure 7). Specifically, our data support the 
notion that GATA4/GATA6 promotes crypt cell prolifera- 
tion by directly repressing Spdefgene expression. 

GATA6 has been suggested to be oncogenic in multiple 
cancers and pre-cancerous lesions, including Barretts 
esophagus, gastric cancer, pancreatic cancer, and colon 
cancer [20]. Gata6 expression is up-regulated in colon 
cancer epithelial cells [21-23], as well as in non-malignant 
cells along the stromal margins in human colorectal cancer 
[22], but it should be noted that it has not been determined 
whether this is a correlative, causative or protective event. 
Recently, Noah et al. [24] showed that SPDEF functions as 
a colorectal tumor suppressor. In colorectal tumors from 
patients, loss of SPDEF was observed in approximately 85% 
of tumors and correlated with progression from normal 
tissue, to adenoma, to adenocarcinoma. Further, SPDEF 
inhibited the expression of p-catenin-target genes in mouse 
colon tumors, and interacted with |3-catenin to block its 
transcriptional activity in colorectal cancer cell lines, 
resulting in lower levels of cyclin Dl and c-MYC. Pathway 
analysis showed that Gata6 deletion also modulated 



c-MYC targets (Additional file 1: Figure SI) and that the 
GATA6/SPDEF overlapping gene set likely functions in 
Wnt signaling (Figure 1). Thus, our data here suggest 
that the possible tumorigenic effects of GATA6 could 
be mediated by its repression of Spdefgene expression. 

Methods 

Analysis of publicly available gene profiling data 

From the National Center for Biotechnology Information 
(NCBI) website (http://www.ncbi.nlm.nih.gov/), we down- 
loaded gene profiling data from the ileum of conditional 
Gata6 knockout mice (GSE22416) [4], and from the small 
intestine of Spdef-/- mice (GSE14892) [9]. Data analysis 
was performed using Cistrome (www.cistrome.org), a flex- 
ible bioinformatics workbench with an analysis platform 
for ChlP-seq and gene expression microarray analysis [25]. 
Significant over-representation of overlapping target genes 
was determined by the Fishers Exact Test. Gene ontology 
analysis and pathway analysis was conducted using Data- 
base for Annotation, Visualization and Integrated Discovery 
(DAVID) [26]. Gene set enrichment analysis was performed 
by using publicly available Gene Set Enrichment Analysis 
(GSEA) tools v2.3 from the Broad Institute [27]. 

Cell culture and lentiviral infection 

Lentiviral infection was conducted as previously de- 
scribed [28]. 293 T and Caco-2 cells were cultured in 
Dulbecco s modified Eagle s medium (DMEM; Mediatech, 
Inc., Manassas, VA) containing glucose (4.5 g/liter), 
1-glutamine, and sodium pyruvate. The media was supple- 
mented with 10% fetal bovine serum (Atlanta Biologicals, 
Lawrenceville, GA) and penicillin (100 U/ml) -streptomycin 
(1 mg/ml) (Sigma Chemical Company, St. Louis, MO). 
Both cell lines were maintained in 5% CO2 at 95% relative 
humidity and 37°C. Media were replaced every 2 to 3 days. 
To produce lentivirus for the knockdown of Gata6, 293 T 
cells were plated at 10^ cells/ml on 6-cm plates and trans- 
fected after 24 h with 1 (ig of vesicular stomatitis virus 
glycoprotein G (VSV-G) envelope-expressing plasmid pMD. 
G, 1 (ig of pCMV-dR8.91 (Delta 8.9) plasmid containing 
gag, poU and rev genes, and 1 (ig of shRNA-pLKO.l 
plasmids expressing a knockdown short hairpin RNA 
(shRNA) (Sigma) for green fluorescent protein (GFP) 
(used as a control [Ctl]) (SHC005) or human Gata6 
{G6kd) using 6 \A of FuGENE 6 reagent (Roche Diagnostic 
Corporation, Indianapolis, IN). A total of five human 
Gata6 knockdown constructs were screened, and only one 
(TRCN0000005392) was subsequently found to reduce 
Gata6 mRNA significantly. Media from the transfected 
293 T cells were changed 16 h after transfection. Culture 
media containing lentiviral particles were collected and 
filtered (0.45 \im pore size) 48 h after transfection and 
transferred immediately to Caco-2 cells plated at 30% con- 
fluence to infect for 2 h using 0.5 \A of polybrene/ml. 
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Infection was repeated the next day followed by selection of 
infected cells using 4 (ig/ml puromycin. Infected cells were 
kept under conditions of selection until the day of har- 
vesting. Trypsinized cells were homogenized using a QIA 
shredder (Qiagen, Inc., Valencia, CA), and RNA was iso- 
lated as described above for mouse tissue. Nuclear extracts 
were isolated as described previously [11]. Gata6 knock- 
down was determined by quantitative reverse transcriptase- 
(qRT-) PGR for human Gata6 (see Additional file 1: 
Figure S3A for primers) and by Western blot analysis as 
previously described [5] using rabbit anti-GATA6 (Cell 
Signaling Technology, Inc, Danvers, MA; Cat. No. 4253S). 

Chromosomal immunoprecipitation (ChIP) assays 

Using a publicly available GATA6 chromatin immunopre- 
cipitation-high throughput sequencing (ChlP-seq) database 
in Caco-2 cells [12], sequences were mapped to reference 
genome Homo Sapiens build 18 (HG18) using ELAND 
tools, allowing 0 to 2 mismatches (Illumina), and binding 
peaks were identified by model-based analysis of ChlP-seq 
(MACS) [29] using default parameters and P value cutoffs 
of 10-^^. 

Chromosomal immunoprecipitation (ChIP) assays were 
conducted as previously described [28]. Caco-2 cells were 
incubated in DMEM containing 1% formaldehyde (Fisher 
Scientific, Pittsburgh, PA) for 10 min at 37°C. The cells 
were washed 2 times with PBS, scraped, and resuspended 
in lysis buffer (50 mM Tris-Cl [pH 8.1], 10 mM EDTA, 
1% sodium dodecyl sulfate [SDS]) containing protease in- 
hibitor cocktail and PMSF (10 (il/ml). The samples were 
sonicated to obtain chromatin fragments of between 400 
and 1,000 bp. Sonicated samples were resuspended in 
ChIP dilution buffer (1% Triton X-100, 2 mM EDTA, 
150 mM NaCl, 20 mM Tris-HCl, pH 8.1) and incubated 
overnight at 4°C with Dynabeads beads (Life Technologies, 
Grand Island, NY) conjugated with rabbit anti-GATA6 
(Cell Signaling; Cat. No. 4253) or rabbit IgG (negative 
control) (Millipore, Temecula, CA; Cat. No. PP64B). The 
IP samples were washed 6 times with radioimmunopreci- 
pitation (RIPA) buffer (50 mM HEPES [pH 7.6], 0.5 M 
LiCl, 1 mM EDTA, 1% Nonidet P-40, 0.7% sodium 
deoxycholate), and the DNA was recovered by reverse 
cross-linking in 1% SDS-0.1 M NaHCOs for 7 h at 65°C. 
DNA was purified using a QIAquick PGR purification kit 
(Qiagen) and quantified by Picogreen (Life Technologies). 
One nanogram of DNA was used per qPGR reaction using 
primers shown in Additional file 1: Figure S3B. 

Mice 

Previously established and confirmed GataS^''''^^^''''^ [30], 
Spdef-/- [9], and transgenic VillinCreER^^ [31] mice were 
used in this study to produce four groups of mice, including 
controls (C^/), conditional Gata6 deletion in the intes- 
tinal epithelium {Gata6AIE), germline Spdef knockout 



(SpdefKO), and double Gata6/Spdef knockout (DKO), 
The genotypes of each group are as follows: 



Ctl: Gata6' 



■loxP/loxP 



, Spdef -\-/-, VillinCreER -negative 



GataGAIE: Ga^^6^^^^^^^^^, Spdef +/- VillinCreER^^-positive 
SpdefKO: GataG^'"'''^^'''''^ Spdef-/- VillinCreER^^-negative 



DKO: GataG 



■loxP/loxP 



, Spdef -/-, VillinCreER -positive 



DNA was extracted from tail biopsies, and genotypes 
were determined by semiquantitive polymerase chain reac- 
tion (PGR) using previously validated primers (Additional 
file 1: Figure S3C). Male and female mice four weeks of age 
were treated with a single dose of tamoxifen (Sigma) 
(0.1 ml, 10 mg/ml; dissolved in ethanol/sunflower oil =1:9 
[vol/vol]) daily for five consecutive days as described [5], 
followed by a single dose two weeks later (see timeline. 
Additional file 1: Figure S4). Mice were killed and tissue 
was collected 28 days after the start of tamoxifen treat- 
ment. Bromodeoxyuridine (BrdU, 0.1 ml of 10 mg/ml) 
was injected two hours before dissection. Approval was 
obtained from the Institutional Animal Care and Use 
Committee. 

Tissue isolation and processing 

Mice were anesthetized for dissection as previously de- 
scribed [28]. The most distal 1.0 cm segment of small 
intestine adjacent to the ileocecal valve was snap frozen 
for RNA isolation, and the 6 cm segment proximal to that 
was removed and transferred to a glass plate on wet ice 
and prepared for sectioning. This segment was flushed 
with ice cold 4% paraformaldehyde (PFA) in PBS, cut 
longitudinally, pinned open onto paraffin wax in a petri 
dish filled with 4% PFA in PBS with the epithelium facing 
upward. After a 5 min incubation period, the pins were 
removed and the proximal end was grasped with forceps, 
rolled with the epithelium facing outwards, and a pin was 
inserted transversely to secure the roll in place. The roll 
was placed in 10 ml of 4% PFA in PBS and mixed gently 
on a tube rotator for 16-18 hr at 4°C. The PFA was 
decanted and the tissue was washed in PBS 3 times for 
20 min each at 4°C, and then dehydrated in 70% ethanol. 
Tissue was processed at the Rodent Histopathology Core 
at the Dana Farber/Harvard Cancer Center (Boston, MA) 
for paraffin embedding and sectioning. Selected slides 
were stained using the periodic acid Schiff (PAS) reaction. 

RNA isolation and gene expression analysis 

RNA was isolated using the RNeasy kit (Qiagen), and 
mRNA abundances were determined by qRT-PCR as 
previously described [4], using validated primer pairs 
(Additional file 1: Figure S3 A). Glyceraldehyde-3-phosphate 
dehydrogenase (GAPDH) mRNA abundance was measured 
for each sample and used to normalize the data. Data were 
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p53, and (B) a decrease in targets of the proto-oncogene c-MYC. (C) 
Legend defining symbols used in the networks. Arrows indicate the direction 
of the interaction. Red circles = up-regulated transcripts; Blue circles = down- 
regulated transcripts. Differentially expressed transcripts were determined at 
the 5% FDR level using Significance Analysis of Microarrays (SAM) and 
interaction networks were developed from the differentially expressed 
transcripts using Metacore. Figure S2. Gata6 and Spdef mRNA abundance 
in ileum in each group of mice. Gata6AIE and DKO mice had significantly 
lower Gata6 mRNA abundances than Ctl and SpdefKO mice. SpdefKO and 
DKO mice had significantly lower Spc/ef mRNA abundances than Ctl and 
Gata6AIE mice. Figure S3. Primers used for: (A) qRT-PCR, (B) ChIP assays, and 
(C) genotyping. Figure S4. Timeline for study. Mice 4 wks of age were given 
Tamoxifen as indicated {black circle) beginning on Day 0, and an 
injection of BrdU 2 hr before tissue collection at Day 28. 



expressed relative to the median value of control ileum. A 
minimum of five mice in each group was analyzed. 

Immunohistochemistry 

Tissue sections were immunostained as previously de- 
scribed [28]. Primary antibodies included rabbit anti-Ki67 
(Thermo Fisher Scientific, Inc., Fremont, CA; Cat. No. RM- 
9106-Sl) (1:200), mouse anti-BrdU (Thermo; Cat, No, 
MS-1058-PO) (1:250), goat anti-cryptdin related sequence 
4C (CRS4C) (gift from Dr, A, J, Ouellette, University of 
Southern California, Los Angeles, CA) [32] (1:2000), and 
rabbit anti-MUC2 (Santa Cruz Biotechnology, Inc, Santa 
Cruz, CA; Cat, No, scl5334) (1:100), Secondary antibodies 
included biotinylated donkey anti-rabbit IgG, donkey 
anti-goat IgG, and donkey anti-mouse IgG (all from 
Vector Labs, Burlingame, CA), Biotinylated antibodies were 
linked to avidin-horseradish peroxidase conjugates (Vector 
Labs), visualized using 3,3'-diamino benzidine (Sigma) for 
2 to 5 min, and lightly counterstained with hematoxylin. 

Cell counting 

The total number of Ki67-, BrdU-, or CRS4C-positive cells 
in crypts was determined as the total number per crypt. 
Only well oriented crypts with the epithelial layer on at 
least one side continuous with the villus epithelial layer 
were counted, and a minimum of 6 crypts per slide were 
analyzed. The average number of CHGA-positive cells 
was expressed as a fraction of total epithelial cells (villi 
and crypts) from a minimum of 5000 epithelial cells per 
slide, with equal representation of crypts and villi. All de- 
terminations were blinded and conducted on a minimum 
of 5 animals per group. 

Statistical analyses 

In Caco-2 cells, qRT-PCR data was compared by the 
students t-test and ChIP data were compared by the ana- 
lysis of variance (ANOVA) followed by the Tukey-Kramer 
multiple comparison test. In mice, mRNA measurements 
had unequal variances across groups requiring nonparamet- 
ric statistics, and were thus compared by the Kruskal-Wallis 
test followed by the Dunn multiple comparison test, and 
presented as individual data points and medians. Cell count 
determinations had equal variances across groups allowing 
parametric statistics, and were thus compared by ANOVA 
followed by the Tukey-Kramer multiple comparison test, 
and presented as mean ± SD, Differences were considered 
statistically significant at a P-value of less than 0,05, 

Additional file 



Additional file 1: Figure SI. Intestinal Gata6 deletion alters gene 
networks controlling cell proliferation in the mature ileum. Network 
analyses on microarray data of Ctl and Gata6AIE ileum (n=3 in each 
group) revealed (A) an increase in targets of the tumor surpressor gene 
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